function P=s_spiral(G,P,k)
P.ang=P.ang+0.5;
P.pos=G.step*[cos(P.ang),sin(P.ang)];
